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ABSTRACT 

High resolution spectra of the Spitzer Space Telescope show vibration-rotation 
absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of 
deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense 
{n > 10^ cm"^), warm (Tex = 200- 700 K) molecular gas with high column densities 
of these molecules ranging from a few 10'^ - 10^^ cm"^. Abundances relative to H2, 
inferred from the silicate optical depth, range from ~ 10"^ to 10"^ and show no cor- 
relation with temperature. Theoretical studies show that the high abundances of both 
C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid 
surrounding an AGN as the origin of this dense warm molecular gas. Galactic mas- 
sive protostars in the so-called Hot Core phase have similar physical characteristics 
with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, 
the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much 
higher toward the (U)LIRGs in the cooler (Tex < 400 K) phase. We suggest that the 
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warm dense molecular gas revealed by the mid-IR absorption lines is associated with 
a phase of deeply embedded star formation where the extreme pressures and densities 
of the nuclear starburst environment have inhibited the expansion of Hll regions and 
the global disruption of the star forming molecular cloud cores, and 'trapped' the star 
formation process in an 'extended' Hot Core phase. 

Subject headings: infrared: ISM - ISM: evolution - ISM: galaxies - ISM: molecules 
- galaxies: nuclei 



Introduction 



One of the holy grails in the study of luminous and ultraluminous infrared galaxies (LIRGs 
and ULIRGs) is to elucidate the true nature of the central energy source. (U)LIRGs emit nearly all 
their energy in the mid and far-infrared part of the spectrum. LIRGs have a luminosity Lg.iooo^m > 
IO^^Lq and ULIRGs L^-mofim > IO^^Lq, equal to t he power output of quasars. (U)LIRGs are gen- 



erally f ound in interacting a nd merging systems (e.g Armus. Heckman. & Mileyll987 : Sanders et al 



1988al : iMurphy et al.lll996l). During the merger large amounts of gas and dust are concentrated 
toward the nucleus (e.g. lMihos & Hemquistill996h . fueling a massive starburst and possibly a mas- 
sive black hole (an AGN). (U)LIRGs are without doubt the most spectacular sites of star formation 
in the Universe and if a substantial part of their energy originates from AGN activity, it would show 
the AGN in its very earliest deeply enshrouded phase. 



The s ensitive Infrared S pectrograph (IRS) ([Houck et al.l l2004) on board the Spitzer Space 
Telescope (iWemer et al.1 120041) has revealed the richness, co mplexity and diversity of t he mid- 
infrare d spect ra toward a sample of deeply obscured (U)LIRGs ( Armus et al.ll2004 , 2006al lbl: rSpoon et al 
2006alJbf) : see iGenzel & Cesarsky daOOQ) for a review of results from the Infrared Space Obser- 
vatory. The general characteristic of these spectra is the presence of deep, broad absorption of 
amorphous silicates, centered at 10 and 18 //m. In addition, the spectra show a large variety in 
absorption features of crystal l ine sili cates, aliphatic hydrocarbons, water ice, and gas phase bands 
of hot CO (ISpoon et al.ll2004ll2006cl) . R^H emission bands are generally weak and in some cases 
absent. Absorption bands of more volatile ices (e.g. CO or CO2), commonly detected in Galactic 
dense molecular clouds, are generally absent or very weak. 

The very compact nuclei of deeply obscured (U)LIRGs are packed with gas. Molecular gas 
has been observed at mil limeter wavelengths — through its low lying pure rotat ional t r ansitio ns — 



in (U)LIRG nuclei (e.g. ISolomon et al.lll997l : lOownes & Solomonlll998h . see I Aaltol (120051) for a 
recent overview. In recent years millimeter emission lines of HCN and HCO "^ have been observed 
revealing the presence of relatively dense (hh ~ 10^ cm"^) molecular gas (|Imanishi et al.l l2006l : 
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Kohnoll2005|) . However, the analysis of the millimeter emission lines is complicated as a result of 
beam dilution which strongly depends on the molecules and transitions observed. 

The IRS Short-High spectrometer (IRS-SH) on board the Spitzer Space Telescope allows for 
the first time the direct study of the very dense and warm molecular gas in the central regions 
of deeply obscured luminous and ultraluminous infrared galaxies through IR pencil-beam line-of- 
sight absorption spectroscopy against the continuum of individual nuclei or unresolve(il3 double 
nuclei. In particular, the mid-infrared vibration-rotation bands of C2H2 and HCN uniquely trace 
warm (100 <T < 1000 K) and dense (nu > 10' cm"^) molecu lar gas. These bands have prev iously 
been detected primarily towa rd Galactic massive protostars (ILacy et al.lll989l : Evans et al.lll991. : 
Lahuis & van Dishoecdbood) . They may prove to be a new strong tool to probe the heating sources 
of deeply obscured (U)LIRG nuclei (starburst or AGN activity). 



2. Observations 



The observations presented in this paper are part of the Spitzer IRS observing programs 
(PID) 105 (IRS GTO ULIRG program, J.R.Houck PI), 96 (IRS GTO program on nearby AGN, 
J.R. Houck PI) and 1096 (Director's Discretionary Time (DDT) proposal, H.W.W. Spoon PI). The 
IRS GTO ULIRG sample comprises ~1 00 ULIRGs in the re dshift range 0.02 < z < 0.93, selected 



primarily from the IRAS 2-Jy sample ([Strauss et al.l 1 19921). IRAS l-J y sample (|Kim & Sanders 



1998h and the FIRST/IRAS radio-far-IR sample ([Stanford et all 120001) . The samples of PIDs 96 
and 1096 contain 3 additional ULIRGs (IRAS 04384-4848, IRAS 03000-2719 and IRAS 021 13- 
2937) and 3 additional LIRGs (IRAS 02530-1-0211, IC 860 and NGC4418). For all sources in 
this combined sample low-resolution spectra (R = A/AA ~ 100) have been obtained, while high 
resolution spectra (R = 600) have been taken only for the brighter half of the sample. 

All high resolution spectra in this sample have been investigated for the presence of vibration- 
rotation absorption bands of C2H2 (13.7 /xm), HCN (14.02 /im) and CO2 (15.0 //m) against the 
nuclear continuum. Fifteen of the sources listed in Table[T]show absorption due to (some of) these 
species. Because of the low signal to noise ratio of the spectra and/or the low H-column densities 
(as derived from the 9.7 um silicate optical depth; see Section 13.31 ). the derived upper limits on 
the molecular column densities for all other sources do not place very stringent constraints on the 
C2H2 and HCN abundances. To illustrate this, four sources without positive detections with both 
a moderate-to-large hydrogen column density and good signal to noise spectra, are included in 
Tables [Hand 121 The upper limits on the C2H2 and HCN abundances fall within the range of derived 



"The IRS-SH slit width is 4.7", equal to the size of the PSF at 19.5^m. 
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abundances toward the fifteen sources. 



Data reduction started from crosstalk corrected echelle images using S12 and S13 Spitzer 
archive data. Pro cessing was done using the c2d analysis pipeline (|Kessler-Silacci et al.l l2006l : 
Lahuis et al.ll2006h . It includes echelle image analysis (a.o. bad-pi xel detection, image arith mic, 
and optimal spectral extraction), defringing of the extracted spectra (|Lahuis & Boogertil2003l) . and 
spectral analysis (see Section [ST]). The optimal spectral extraction uses an analytical source profile 
defined and calibrated using a suite of calibrat or stars. Calibrati on is done using MARCS stellar 
models provided by the Spitzer Science Center (|Decin et al.ll2004t) . The optimal spectral extraction 
employs a combined source and local sky fit. This provides wavelength dependent sky estimates 
and allows to discriminate between resolved and unresolved spectral features. 

Figure [H shows the continuum divided, sky and redshift-corrected spectra from 13.5-15.5 /im 
covering the absorption bands of C2H2, HCN and CO2. Plotted in red are the best-fit single temper- 
ature synthetic spectra (see Section [3]). Indicated with the dotted verticals are the positions of the 
ground-state Q-branch transitions of the three molecules. For IRAS 15250-1-3609, IRAS 20100- 
4156, and IC 860 a small section of the spectrum in between HCN and CO2 is affected by artifacts 
in inter-order sections of the spectrum, and these sections have been clipped from the presented 
spectra. 



3. Analysis 
3.1. Molecular analysis 



The Spitzer spectra unambiguously reveal the presence of the Q-branch transitions of C2H2 
z/5 = 1 - 0, HCN z/2 = 1 - and CO2 1^2 = 1-0, each of which blends into a "broad" absorption 
feature. The corresponding P- and i? -branch transitions of these species are difficult to observe with 
the Spitzer-\RS) due to spectral dilution at the IRS resolving power and the presence of instrumental 
fringe residuals with frequencies and amplitudes close to those of the P- and 7? -branch lines. The 
2-branch transitions are analyzed using a pure absorption model assuming local thermodynamic 
equilibriur n (LTE) excitation of the levels a t a si ngle temperature. The adopted method is described 
in detail in iLahuis & van Dishoeckl (120001) and lBoonman et al.l (120031) . which includes references 
to the molecular parameters and data used in the model. The main fit parameters are the excitation 
temperature and column density along the line of sight for a given intrinsic line width, defined by 
the Doppler Z^-value. It is assumed that the absorbing molecules have a covering factor of unity 
of the continuum, i.e. the mid-IR continuum is composed solely by the regions toward which 
the molecular absorption features arise. A covering factor less than unity, i.e. larger continuum 
emitting regions, increases the true optical depth of the absorption features resulting in higher 
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column densities and possibly lower temperatures (when lowest excitation lines saturate at large 
column densities). 

The derived excitation parameters do not strongly depend on the exact value of b. Only for 
low values of b (< 5kms"') will saturation result in an underestimate of the column density for 
the sources with the largest column densities. Such low b values are expected for quiescent gas 
where thermal broadening dominates. However, non-thermal broadening will likely dominate in 
the dense energetic ISM of the galactic nuclei and larger b values are expected. 

A direct estimate of b is obtained from spectr ally resolved CO u =1-0 absorption lines 
toward the north-west nucleus of IRAS 08572-1-3915 (|Geballe et al.ll2006|) . which shows a complex 
velocity structure. A low column density, cold CO component absorbs near the systemic velocity. 
The CO absorption is however dominated by a broad (FWHM ~ 200kms~^) blue-shifted warm (> 
200 K) gas component. This is most likely the conglomerate of a velocity distribution over multiple 
lines of sight within the beam. The continuum confusion then requires a higher column de nsity 
than that of 2 x 10'^ cm"^ estimated from the observed optical depth by lGeballe et al.l (120061) . We 
have fitted synthetic LTE absorption spectra to the higher excitation CO lines using the spectrum of 
Geballe et al.l (|2006h but with the spectral resolution reduced to match the observed profile width 
{R ~ 1500). This allows to estimate a velocity averaged column density for a given b value. Values 
< lOkms"' result in progressively poorer fits. Good fits are made for b values of 10-25kms~^ 
requiring column densities of ~ 2 x 10'^ cm"^ down to ~ 5 x 10'^ cm"^. The Doppler b-va\ue has 
therefore been fixed to 20km s~' for all sources. 



Recent observations using TEXES (|Lacy et al.ll2002|) . visitor instrument on Gemini-North, 
have revealed spectrally resolved (FWHM ~ 80km s~^ ) blue-shifted absorption in the C2H2 1^5 = 1-0 
R(13) ro-vibrational line toward one of our sources, NGC4418 (Knez, private communication). 
The current data do not allow to put more stringent constraints on the Doppler Z7-value. However, 
the spectrally resolved C2H2 absorption could be an indication that multiple (spatial and velocity) 
components of the warm and dense gas are common in (U)LIRG nuclei. 



3.2. Fit results 

Table |2] lists the derived excitation temperatures and column densities from best fit synthetic 
spectra to the continuum divided (U)LIRG spectra. Excitation temperatures ranging from 200 to 
700 K and column densities of a few x 10'^ to 10'^ cm"^ are observed. These results are derived 
from a simultaneous fit to the three absorption bands, in which the excitation temperature, set to 
be the same for all three molecules, is constrained by the Q-branch profiles of C2H2 and HCN. 
Since the analysis is restricted to that of the analysis of the resolved Q-branches (see 13.11) . the 
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derived excitation temperatures are not well constrained (~ 30%). However it does clearly allow to 
discriminate between molecular gas with warm (~ 200-300 K) and hot (~ 500-700 K) excitation 
temperatures. This is illustrated in Figure [21 by fits to the observed absorption profiles of IC 860 
and IRAS 01003-2238, for which we derive excitation temperatures of 280 and 630 K respectively. 

The 2-branch profiles of C2H2 and HCN reveal significant contributions to the absorption 
from relatively high rotational levels. Hence, given the high critical density of these levels, the 
absorbing gas is likely very dense (n > 10^ cm"^). Observations of the ro-vibrational transitions of 
CO in IRAS 08572+3915 and IRAS 00 183-71 1 1 indeed imply densities in excess of 3 x 10^ cm'^ 
dOeballe et allbood : Ispoon et al.l[2004h . 



3.3. Abundances 



The derived column densities are translated into abundances by using a total hydrogen column 
density obtained from the apparent optical depth of the 9.7yum silicate absorption band listed in 
Tabled! The optical depth is converted assuming: 

A^H = r9.7(3.5 X 10^^)cm-2, 



appropriate for the dust in the solar neighbourhood (|Roche & Aitkenll 198411 1985b . 



This relation is most appropriate for embedded (extinction dominated) sources probing the 
absorption along the pencil beam line of sight toward the warm central continuum emitting region. 
To illustrate, for the sample of Galactic massive protostars (see Section |431) the hydrogen column 
density derived using the silicate optical depth estimated from the ISO-SWS archive spectra agrees 
within a factor of two with the hydrogen column de nsity derived from CO = 1 - measurements 
dMitchell et al.lll990l : Lahuis & van Dishoeckll200ol) . 



Another estimate of the applicability of this method is provided through the analysis of the 
4.6 fxm CO u = 1-0 absorption spectrum toward IRAS 08572-1-39 15NW, which gives CO column 
densities of ~ 5 x 10'^ cm"^ up to ~ 2 x 10'^ cm"^ depending on the assumed value of b (see 
Section 1311). Adopting Nqo/Nu = 10~^ (all gas-phase carbon in CO) this resuks in A'^h = 5 x 10^^- 
2 X 10^^ cm"^. This is in reasonable agreement with the hydrogen column density derived from 
the silicate optical depth (A^h ~ 1-5 x lO^^cm'^). If anything, the IRAS 08572-1-39 15NW data 
suggest the hydrogen column density might be slightly overestimated (e.g. as a result of foreground 
extinction) resulting in underestimating the derived molecular abundances. We assume a similar 
uncertainty for the other sources. 

The derived abundances and H2 column densities (assuming = Nu/2) are listed in Table 
|2l The HCN abundance is up to three orders of magnitude higher than in cold molecular clouds 
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(IBergin & Langerlll997|) . High abundances of HCN, C2H2, and CO2 have also been measured for 
Galactic massive protostars (see Section |431 ). Figures[3]and|4]compare the abundances to chemical 
models and results from studies of Galactic massive young stellar objects (Section IH). 



3.4. Gas temperature 

The derived excitation temperatures range from 200 to 700 K. While statistics are small, of the 
seven sources with the highest gas excitation temperatures {T^^x = 400-700 K) five sources show the 
highest 5.5 fim to 28 fim continuum flux ratio. These galaxies also show a rising near-IR continuum 
typical for hot dust while the other galaxies show a decreasing near-IR continuum characteristic 
for emission dominated by stellar photospheres. It seems thus that there is a connection between 
the observed dust and molecular excitation temperatures. 

The presence of absorption (rather than emission) lines implies a strong gradient in the tem- 
perature and an underlying hot-dust continuum against which the absorption lines are formed. 
However, the derived excitation temperatures of the molecular gas for the warmest sources {T^x = 
400-700 K) are higher by a factor of 2 to 3 than the mid-IR color temperature derived from the 5.5 
to 28.0 fim flux ratio. While the discrepancy between the excitation temperatures and the mid-IR 
color temperature may (in part) reflect the uncertainties in the analysis, the observed colder dust 
continuum may point towards the importance of continuum extinction by (a shell of) colder fore- 
ground dust exterior to the C2H2 and HCN absorbing layer, or the predominance of emission of 
cold nuclear dust within the IRS aperture (i.e. a small beam filling factor for the warm dust and 
gas). Both may result in an added uncertainty in the derived abundances resulting from overesti- 
mating the hydrogen column density or from a covering factor of the molecule rich warm gas less 
than unity (see Section [3T| ). 



4. Discussion 



4.1. Warm molecular gas in (U)LIRGs 



The Spitzer-IRS spectra reported here reveal surprisingly strong mid-IR absorption bands of 
C2H2, HCN and CO2 toward a sample of (U)LIRG nuclei. These absorptions reveal the presence 
of c opious amounts of warm molecular material an d are in line with earlier detections of warm CO 
gas (|Spoon et al.ll2004 l2006al : iGeballe et al.ll2006|) . This molecular gas is very dense (n > 3 x 10^ 
cm"^) and warm (r ~ 200-700 K). Given the column density and estimated density, this gas 
occupies only a small fraction of the nuclear region (~ 0.01 pc) and, given its high temperature, is 
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likely located near the intrinsic mid-infrared source of luminosity of these regions. 



At present there are still substantial uncertainties in accurately inferring the nature of the 
power source domin ating the bolometric luminosity in infrared luminous galaxies. As discussed in 
Armus et al.l (l2006bl) . it appears that there are disagreements among various diagnostic diagrams 
based on line ratios, continuum colors or PAH equivalent width. Perusal of the 2-40/im spectra 
reveal that the sources with molecular gas absorptions appear to be a "mixed bag". While most 
sources show a deep silicate feature, some (e.g., Mrk231, IRAS 01003-2238, IRAS 05189-2524 
and IRAS 13218-1-0552) do not. Indeed, the 2-40/im continuum spectra of these sources show 
strong similarities with the spectra of bona fide AGN-dominated spectra with shallow silicate fea- 
tures. In contrast, sources like IRAS 17208-0014 and IRAS 22491-1808 look quite starburst-like. 
Finally, sources such as IRAS 08572-1-3915 show neither signs of starburst activity (e.g., PAH fea- 
tures) nor signs of AGN activity (e.g., high ionization lines or a broad Ha line). While the presence 
of hot dust, as visible in the near-IR, may be taken as evidence for the presence of an embedded 
AGN in this source, any deeply obscured source - including deeply embedded proto stars - will 
produce hot dust. Indeed, the overall mid-IR spectrum of IRAS 08572-1-3915 resembles that of 
many galactic massive protostars, such as AFGL2591, with a preponderance of dust absorption 
features. 

The C2H2, HCN, and CO2 molecular lines may provide a new probe of the conditions in the 
inner deeply obscured nuclei and hence shed light on the ultimate power source of (U)LIRGs; e.g., 
starburst phenomena versus AGN activity. In particular, the presence of copious amounts of warm 
molecular gas enriched in C2H2, HCN, and CO2 is common in regions of massive star fo rmation 
dsoonman et al.ll2003l : Lahuis & van Dishoecklboool : Lacv et al.l[l989l : Isvans et al.lll99lL while, 
on the other hand. X-rays associated with a central black hole can influence molecular chemistry 
and abundances far beyond the mean free path length of stellar (E)UV photons. The molecular 
observations are discussed in the light of these scenarios in the next subsections. While the possible 
errors associated with the derived abundances are recognized, the discussion will center on Figures 
mandU 



4.2. AGN activity and X-ray driven chemistry 



Sub-millimeter studies have revealed high HCN/CO (and HCN/HCO"^) abundance ratios in 
Seyfert type I nuclei and this has been taken as an indication of t he importance of X-ray illu- 
minat ion for the composition of the molecular gas in these nuclei (|Imanishi et al.l l2006l : iKohno 
20051). Some of the sources in our sample are known to harbor an AGN. In particular, Mrk 231 and 
IRAS 13218-1-0552 are QSOs, IRAS 08572-1-3915 and UGC5101 are LINERs, and IRAS 01003- 
2238, IRAS 05189-2524 and Mrk 273 have optical LINER/Seyfert 2 spectra (see Table©. Hence, 
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the observed high abundance of HCN, and by inference C2H2 (which has no pure rotational transi- 
tions), may indicate the presence of a buried AGN in all of our sources. The complex geometry and 
inhomogeneous dust distribution associated with the central toroid may then preclude our detec- 
tion of the putative X-rays photons from the central engine driving the chemistry of the molecular 
gas. The mid-IR observations presented in this paper add to the submillimeter results the predom- 
inance of warm molecular gas. Within the AGN heating scenario, it seems obvious to attribute 
the high temperature of the molecular gas to the importance of X-ray heating in a so-c alled XDR 
(IMaloney et al.lll996l : lMeijerink & Spaansll2005|) . Alternatively iGao & SolomonI (|20o4 argue that 
most (U)LIRGs are dominated by starburst power rather than AGN power based on the HCN 
7=1-0 luminosity and CO 7=1-0 luminosity relation. In a starburst scenario the HCN/CO ratio 
traces the amount of dense molecular gas in the galactic ISM. Our observations may present an 
additional test for these two scenarios. 

X-rays are not only an efficient heating agent of the gas but also have a major influence on th e 
chemical composition. Thi s has been explored in the context of AGNs bylMalonev et al.l (119961) : 



Meijerink & SpaansI (|2005|) . and for regions around pro to stars by iStauber et al.l (|2005l . 120061) 7 At 
the high densities relevant for these regions, the gas temperature is set by a balance between X-ray 
heating and gas-dust coUisional cooling. As an example, a gas temperature of 200 K at a density 
of ~ 3 X 10^ cm"^ requires a X-ray flux of ~ 30 erg cm"^ s"^ which for a typical Seyfert galaxy 
with an X-ray luminosity of 10"^"^ erg s"' corresponds to a distance of 50 pc for an unobscured view 
of the central engine. For the molecular features to be seen in absorption, the mid-IR continuum 
has to arise from an inner, warmer zone (e.g., exposed to a higher X-ray flux). Since most of 
these (U)LIRGs show no evidence for strong X-ray emission while the column density associated 
with the molecular gas (10^^- 10^^ cm"^) is small compared to X-ray attenuation lengths, in this 
AGN-scenario, the direct view of the nucleus would have to be blocked by a high-column-density, 
edge-on toroid while the warm molecular gas is offset to the polar regions and has a little-obstructed 
view of the nucleus. 

An AGN origin for the high abundance of warm HCN, C2H2, and CO2 in these (U)LIRGs 
however faces severe theoretical difficulties. Specifically, the high X-ray flux requir ed to explain 
the observed temperatures readily dissociates the molecules and drives the gas atomic ( Maloney et al 
19961 : iMeijerink & Spaansll2005[) . Indeed, calculated HCN abundances are typically less than 10"^ 
wrt. H in X-ray illuminated gas while C2H2 is virtually destroyed for X-ray fluxes in excess of 1 erg 
cm"^ s"^. This is illustrated in Figure [3] which presents the static steady-state chemical composition 
after equilibrium is reached (~ 10^ yr) for clouds illuminated by X-rays. Here, static means that 
the chemistry has been evolved in time - while the X-ray flux, tem perature, and density are fi xed - 



ipera 

) 1stai 



until steady state is achieved in the chemical abundances (see also IStauber et al.ll2005L 120061) . The 
X-ray spectrum is assumed to be thermal in nature with 7x=3 x 10^ K. The ch emistry is relatively 
insensitive to the assumed shape of the spectrum (e.g., IMaloney et al.lll996l) . As these models 
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demonstrate, strong X-ray irradiation leads to decreased abundances o f molecular species such 
as HCN and C2H2. These results are consistent with earlier studies (e.g. ILepp & Dalgarndll996() 
which show an initial increase in the HCN abundance with increasing ionization rate / nu up to 
10"^*^ cm"^ s"^) but then an effectively complete destruction of HCN for (^/^h > 10"^' cm"^ s"^ For 
comparison, the X-ray flux implied by a gas temperature of 200 K (~ 30 erg cm"^ s"') corresponds 
to an ionization rate of 5 x 10"^' cm"^ s"^ at a density of lO** cm"^. Clearly, the warm molecular 
gas revealed by these observations is inconsistent with strong X-ray illumination. 



4.3. Static hot-core chemistry 



At first sight our mid-infrared results resemble those toward a sample of massive Galactic 
young stellar objects (YSOs). For these, the molecular absorption features originate in the warm . 



dense gas near the newly formed star during the so-called hot cor e phase (IBoonman et al.l 12003 



Lahuis & van Dishoecld I2OOOI : iLacy et al.lll989l : Evans et al.lll99lh . This is generally thought to 
represent a very early phase (At <10^ yr) in the formation of massive stars, during the transition 
from the deeply embedded stage to the phase where ionizing photons can escape the protostellar 
accretion envel ope and create first a hyper- and later an ultra-compact Hll region (see review by 
Cesaronil |2005|) . The molecular composition of the hot core is very different from that of cold 
molecular clouds. This is thought to reflect the evaporation of ices - resulting from accretion 
and surface chemistry duri ng the cold preceding dark cloud phase - when the envelope is heated 
by the newly formed star (^ Walmsley & Schilke 1993). Subsequent high temperature gas-phase 
chemistry significantly enhances the abundance of e.g. CH4, C2H2 and HCN up to three orders 
of magnitude, compared to cold dense molecular cloud abundances , for the most evolved sources 
Jviti&WilliamJ 19991 : lootv et al.ll2002l : iRodgers & Chamlevll2003l) . 



In Figure|4]the observed abundances and abundance ratios toward the (U)LIRGs are compared 
to those of a sample of Galactic massive YSOs. Both the YSO and the (U)LIRG sample show 
comparable, large variations in gas temperature and abundances. However, the HCN and C2H2 
abundances show a positive correlation with the gas temperature for the YSO sample but not for 
the (U)LIRG sample. Indeed, at low temperatures (< 400 K), the (U)LIRG abundances are some 
two orders of magnitude larger than those in YSO spectra. 

For the Galactic YSOs, the observed correlation between the excitation temperature and the 
observed abundance s (Figure S]) are in fair agreement with chemical models for such regions 
(|Stauber et al.ll2005|) . The absence of a similar temperature correlation for (U)LIRGs may reflect 
a systematic (and large) error in the abundance determination (Section [33l) or it may reflect a true 
difference between the chemical or physical evolution of the warm gas of the YSOs and (U)LIRGs. 
In this respect, we note that the C2H2/HCN ratio for both samples are very similar, however the 
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CO2/HCN and CO2/C2H2 ratios are reduced for the (U)LIRG sample up to an order of magnitude. 
If static steady-state hot-core chemistry were to apply to the (U)LIRGs as it does to the YSOs, CO2 
would have been detected toward most of the sources in our sample. The fact it has not, suggests 
that a different chemistry may apply to the warm molecular gas in these (U)LIRGs. 



4.4. Pressure confined starburst chemistry 

Static steady- state hot-core and X-ray chemical models have difficulties to consistently re- 
produce the observed high abundances of warm C2H2 and HCN gas and low abundance of CO2 
gas for most of our (U)LIRG sample. Essentially, static chemistry has a hard time producing high 
abundance of HCN and C2H2 in the warm (100-400 K) gas characteristic for most of our (U)LIRG 
sample, while at the same time reducing the CO2 abundance. One possible solution to this conun- 
drum is to produce high C2H2 and HCN, and low CO2 abundances in hot gas (> 800 K), and then 
transport it outward to the cooler gas on a dynamical timescale which is rapid compared to the 
chemical timescale. At high temperatures, the hydrocarbon and nitrogen chemistry are enhanced 
as most of the oxygen is converted into water by neutral-neutral reactions. The abu ndances of 



molecules such as C2H2, CH4, an d HCN can be increased by orders of magnitude (e.g. iDoty et al 



20021 : iRodgers & Chamleyll2003|) while at the same time the formation of CO2 is reduced as its 
primary formation route through OH is blocked. For cold Galactic molecular clouds, the chemical 
timescale is set by the cosm ic ray ionization rate and is about 3 x 10^ yr, independent of density 



(|Glassgold & Langeiill973|) . In a (U)LIRG (or AGN) environment, the cosmic ray flux may be 
increased considerably, shortening this timescale. Likewise, neutral-neutral reaction channels may 
open up in warm gas, further reducing the chemical timescale. In any case, these timescales are 
much shorter than the evolutionary timescale of massive stars or of the starburst. With the typical 
sizescale of the warm molecular zone (~ 3 x 10'^ cm), a timescale of 3 x 10^ yr translates into a 
'diffusion' velocity of only 0.03 km s"' . Even a 100 times faster chemical timescale only requires a 
'diffusion' velocity of a few km s"^ Since the chemical models are well capable of explaining high 
abundances for both C2H2 and HCN and a low abundance for CO2 at high temperatures and den- 
sities, and in light of the discussion above, this 'mixing' is attributed to the global activity created 
by pressure confined massive starburst activity. The strong gravitational potential in the nuclei of 
these galaxies inhibits the disruption of the surrounding warm molecular envelopes by Hll regions 
and supernovae, while producing sufficient turbulent (or wind) motion to distribute warm C2H2- 
and HCN-rich and C02-poor gas in the colder outer envelope regions. 



5. Summary 



We have observed the absorption lines due to the Q-branch transitions of C2H2, HCN, and 
CO2 in the mid-IR spectra of a large number of (U)LIRGs. These observations reveal for the first 
time the presence of copious amounts of warm (200-700 K), dense (n > 10^ cm"^) molecular gas 
in these nuclei. The origin of this warm molecular gas is unclear. Theoretical models show that the 
X-ray fluxes implied by the elevated gas temperatures rapidly destroy these molecules and hence 
it is unlikely that this warm molecular gas is associated with a dense toroid surrounding an active 
central engine. Warm, dense gas rich in C2H2, HCN, and CO2 is commonly observed towards 
galactic massive protostars and is associated with a short-lived phase - the Hot Core phase - be- 
fore advancing ionization fronts disrupt the core. The high molecular abundances in galactic Hot 
Cores are well understood chemically. However, the derived abundances of C2H2 and HCN for the 
cooler (U)LIRGs (T^^ ~ 200-400 K) in our sample as well as the C2H2/CO2 and HCN/CO2 ratios 
are very different from those in galactic protostars. We suggest that this warm molecular gas is 
associated with a phase of deeply embedded star formation in (U)LIRGs where the high pressures 
and densities have inhibited the disruption of the star forming, molecular cores by prohibiting the 
expansion of HII regions, trapping the star formation process in an extended Hot-Core phase. The 
chemical differences between these (U)LIRGs and galactic Hot Cores may then reflect enhanced 
mixing between warm and cold phases due to the high degree of turbulence or wind motion asso- 
ciated with extreme starburst environments. Pressure confined massive starburst activity may thus 
be the driving energy source behind the observed C2H2 and HCN rich warm molecular gas and 
responsible for most of the near-IR characteristics of the deeply obscured (U)LIRGs. 

The current analysis is predominately based on the observed abundances derived from moder- 
ate resolution, mid-IR vibration-rotation absorption bands of C2H2 and HCN. This analysis can be 
extended using high resolution, velocity resolved, ground-based studies of the P- and R branches 
of these molecules as well as from the = 1 - rotation-vibration band of CO. In addition, ob- 
servations of more molecular rotation-vibration bands and high excitation submillimeter lines of 
a large sample of molecules may be instrumental for progress in this field. Specifically, velocity 
and spatially resolved infrared and millimeter data can constrain the source morphology. In com- 
bination with optimized physicochemical models for individual sources it may become possible to 
draw firm conclusions about the physical characteristics of the warm molecular gas and the true 
nature of the power sources in the (U)LIRG nuclei. 

The authors would like to thank Tom Geballe for sharing early results and the CO data of 
IRAS 08572-I-3915NW; John Lacy and Claudia Knez for sharing the TEXES data; and Bernhard 
Brandl, Kees DuUemond, Masa Imanishi, David Rupke, and Marco Spaans for many useful dis- 
cussions. Astrochemistry in Leiden is supported by a Spinoza grant from NWO. 
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"Ramp integration time x number of cycles X number of slit positions 

''Opti cal/NIR spectral classifications taken from: (l)'Armus. Heckman. & Milev ('l989|), (2l lSanders etai]h988bh . r3') IVeilleux etalJh995l) . (4) 
iLow et al. ( 1988), (5) Darling & Giovanelli (2002), (6) Veilleux, G oodricli. & HilU1997n 

'assuming Ho = 71 kms"' Mpc"' , = 0.27,f2A =0.73,r2K =0 



Table 2. LTE excitation parameters and spectral characteristics 
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^The excitation temperature is poorly constrained and can be uncertain up to 30%. See Section lT2l and Figure|2]for details. 
''Abundances wrt. H2 assuming Nw = rg 7(3.5 X 10^") cm"" and Nw = 2 X Nw-^ (see Section l33t . 
''Continuum at 14 /xm in the rest wavelength frame. 

'^Signal to Noise Ratio (SNR) estimated from the residuals after subti'action of the synthetic spectrum. It varies over the covered wavelength 
range. 

'Temperature fixed to 300 K for derivation of column density upper limit estimates. 
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Fig. 1 . — Continuum divided Spitzer-IRS spectra of a sample of (U)LIRGs showing the absorption 
bands of C2H2 and HCN and some of CO2. Plotted in red are best-fit synthetic spectra assuming 
a single excitation temperature for all three molecules. All spectra have been shifted to the rest 
wavelengths. 
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Fig. 2. — Illustration of the temperature sensitivity of the C2H2 and HCN Q-branch profiles in di- 
rect comparison to the observed spectra of IC 860 (best fit T^^ = 280 K) and IRAS 01003-2238 (best 
fit Tex = 630 K). The shape of the Q-branch profile determines the derived excitation temperatures 
with an uncertainty of ~ 30%. A better constrained error estimate is difficult to determine at the 
IRS resolution. 
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Fig. 3. — Results of equilibrium X-ray enhanced chemical models for the three observed 
molecules. The calculated abundances are shown for moderate (10^ cm"^) and high density 
(10^ cm"^) and four X-ray flux levels (Fx = 0, 10"^, 1 and 10^ erg cm"^ s"' plotted with solid, 
dashed, dot-dashed, and long-dashed lines). Note that except for HCN and CO2 at high density, 
the abundances for the highest X-ray flux are belo w the plotted ra nge. Qverpl otted are the observed 
abundances. Details of the model can be found in lStauber et al.l ( 20051 12OO6 ). 
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Fig. 4. — Top: observed C2H2, HCN and CO2 abundances using a total H2 column obtained from 
the 9.8/xm silicate absorption band. Bottom: abundance ratios of C2H2, HCN and CO2. Both 
the absolute abundances and the abundance ratios are indicated with circles and are presented as 
functions of excita tion temperature. Inc l uded are the observed values toward a sample of Galac- 
tic massive YSOs ( Boonman et al. 2003 : Lahuis & van Dishoeck 20o d) plotted with grey triangle 
symbols. The derived absolute abundances toward the (U)LIRGs show significant enhancements 
in the abundances of C2H2 and HCN for sources with excitation temperatures below 500 K. The 
C2H2/HCN abundance ratios are quite similar, however the CO2/C2H2 and CO2/HCN ratios toward 
the (U)LIRGs are systematically lower than toward the Galactic massive YSOs suggesting that a 
different chemistry may apply. 



